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Abstract The single bipolaron problem is examined in the context of the ID Holstein-Hubbard model, em- 
phasizing analogies and differences with respect to the complementary single polaron physics. The bipolaron 
band structure below the phonon threshold is revealed, showing a complex relationship between numerous 
excited bands as the adiabatic limit is approached. Light bipolarons with significant binding energy, the 
stability of large bipolarons, the small to large bipolaron crossover as a function of the Hubbard repulsion, 
as well as the bipolaron dissociation, are investigated in detail, disentangling adiabatic, nonadiabatic and 
lattice coarsening effects. It is emphasized that condensation of bipolarons occurs in the dilute limit only 
at very low temperatures. 
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1 Introduction 

The polaron represents a quasi-particle involving a sin- 
gle electron coupled to the lattice degrees of freedom. 
This coupling introduces local correlations between the 
electron and the lattice field, characterized by a finite 
electron-phonon correlation length dpoi and an effective 
mass nipoi- The translational symmetry of the lattice is 
preserved because the electron and the lattice field (de- 
formation) can travel together. When the local Coulomb 
repulsion is not too strong, a pair of electrons can bind 
by sharing a common lattice deformation field. Such an 
entity is called a bipolaron. The bipolaron condensation 
energy Ai,p is defined as the difference in binding energies 
of the bipolaron Ehp and two uncorrelated polarons Epoi , 
Afip = 2Epoi — Ebp ■ The bipolaron mass mi,p is associated 
with a joint motion along the lattice of a pair of electrons 
and the corresponding lattice field. 

The present work is focused on low-frequency proper- 
ties of a single bipolaron in the context of the Holstein- 
Hubbard (HH) model [2] . Previous investigations of bipo- 
larons based on this model range from variational studies 
|3I4I5I6I7I5] to various numerical approaches, such as exact 
diagonalization calculations for finite clusters 19110111112] 
and the infinite lattice [13114) . quantum Monte-Carlo sim- 
ulations [15116117] , the density matrix renormalization group 
[IS] , and the cluster perturbation theory [12] . Generaliza- 
tions of the Holstein-Hubbard model have been studied 
[20121] by varying the range of the electron-phonon inter- 
action. Perturbative expansion around the atomic limit is 
analyzed comprehensively in reference [22] . whereas the 
exact Green's function for the two-sites problem is de- 



rived analytically in reference [23] . A broader overview of 
the bipolaron related physics may be obtained from refer- 
ences El], El] and US] . 

Although a number of methods have been developed in 
the last decade to investigate the HH bipolaron problem, 
the description of an important part of the phase diagram, 
which involves large adiabatic bipolarons and associated 
crossovers to other regimes, is still scant. The problem lies 
in electron-phonon correlations that emerge when the cou- 
pling is strong and the electron-phonon correlation length 
dbp is large, d^p/a ^ 1 (a is the lattice constant). Namely, 
such correlations impose particular demands on an accu- 
rate, quantum treatment of bipolarons since the corre- 
lations with many phonons at large distances from two 
electrons have to be taken into account. 



For exact diagonalization approaches [13114127] that 
are based on a truncation of the Hilbert space, the num- 
ber of states that should be considered grows exponen- 
tially with increasing coupling and increasing dbp- This 
limits the accessibility of the phase diagram for such cal- 
culations. Indeed, to best of our knowledge, all present nu- 
merically accurate studies of the quantum HH bipolarons 
are restricted either to weak couplings or to limited val- 
ues of the adiabatic ratio t/ojQ < 2, where t is the electron 
hopping energy and ojq is the optical phonon energy of 
the HH model. In other words, the formation and prop- 
erties of large bipolarons with significant binding energy 
still represents a challenging subject. 

For the half-filling case, the bipolaronic phase has been 
studied by the dynamical mean field theory (DMFT) [55] . 
However, this particular case of high charge concentra- 
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tion necessarily differs from the dilute limit, when the 
single-bipolaron theory applies. Furthermore, since its dia- 
grammatic expansion neglects vertex corrections involving 
phonons at different lattice sites, the DMFT for finite- 
dimensional systems fails to describe properly the adia- 
batic electron-lattice correlations spanning several lattice 
sites [55]. For the HH model, with local coupling and lo- 
cal phonons, this problem is most pronounced for the one- 
dimensional D ^ 1 system because, in this case, large adi- 
abatic (bi)polarons are stable, while such correlations for 
higher dimensional D > I systems |30l31j are short-lived. 

The treatment developed here combines the results 
found in the adiabatic limit with numerical results ob- 
tained by the recently proposed relevant coherent state 
method (RCSM) [32]. Within the RCSM, (bi)polaron states 
are obtained by solving a generalized eigenstate problem. 
The latter is defined by choosing a trial set of the most 
relevant wave functions after a careful analysis of the prop- 
erties of the low-frequency adiabatic and nonadiabatic cor- 
relations. 

The RCSM offers an improvement over previous cal- 
culations in several different ways. First, it can be ap- 
plied to any value of the adiabatic ratio t/ujQ, including 
the t/ujo ^ 1 part of the bipolaron phase diagram where 
the strongly-coupled large bipolarons emerge. Second, the 
method provides the full low-frequency band structure of 
the bipolaron states including the excited coherent bands, 
which, to the best of our knowledge, have not been pre- 
viously reported. Such analyses allow a detailed descrip- 
tion of the bipolaron low-frequency dynamics and extend 
the previous investigations, which were mainly focused on 
the properties of the ground state, to include the effec- 
tive mass and, in some cases, the dispersion of the lowest 
bipolaron band. The existence of the first excited zero- 
momentum state below the phonon threshold for inelastic 
scattering had been demonstrated in reference |20) . 

The importance of a successful treatment of the long- 
range adiabatic correlations can easily be set in a broader 
context, involving models for which the electron-phonon 
coupling is not purely local (on-site) as in the HH model. 
Namely, for any dimension D, the increasing range of the 
electron-phonon interaction results in an increment of the 
correlation length dtp- For D > 1, this increment necessar- 
ily introduces large adiabatic bipolarons into the phase di- 
agram. In this respect the RCSM, used here for _D = 1, has 
additional advantages since it can be implemented with a 
high accuracy to a broad range of models and D > 1. 

The present paper is organized as follows. After a short 
introduction to the problem in Section [2j Section [3] cov- 
ers the single bipolaron case, starting with the adiabatic 
limit and emphasizing analogies with adiabatic polarons. 
In the next step, differences between the polarons and 
bipolarons appearing with nonadiabatic correlations are 
examined. This is followed by an analysis of the role of 
the Hubbard repulsion U on the condensation energy for 
both large and small adiabatic bipolarons. In Section|3]the 
RCSM is used to calculate accurately the bipolaron spec- 
tra below the phonon threshold for the incoherent scatter- 
ing. Beside the lowest coherent bipolaron band, additional 



coherent excited bands are found to be related to excita- 
tions of adiabatically softened phonon modes of the mov- 
ing lattice field. Particular attention is paid to the small 
and large U adiabatic regimes and the parameter space in 
which small light bipolarons are formed with substantial 
condensation energies. Section [5] gives a brief discussion 
of finite temperature effects, with emphasis on particular 
aspects of the dilute limit. A summary of the results is 
given in Section [51 

2 General 

The ID Holstein-Hubbard Hamiltonian is given by [5] 

n 71, s 

^ 2ff ^ ^ '^n.s'^n,s U„ + U ^ ^ ^n,s'^n,scl^ ,,iC„,s' : (1) 
n,s n,s^s' 

where Cn,s is the annihilation electron operator on site n 
with spin s =t, ]r, Un is the dimensionless lattice displace- 
ment for the site n, and T represents the lattice kinetic 
energy shifted, for convenience, by the zero-point energy 
of the free lattice. 




The Hamiltonian (|T]) describes electrons in the tight bind- 
ing nearest-neighbor approximation coupled to the disper- 
sionless branch of optical phonons. The electron-phonon 
and electron-electron interactions, given by the last two 
terms in equation ([T]), are local. While the latter inter- 
action is instantaneous, the former contains retardation 
effects through the lattice coordinates u„. Un = 1/2 cor- 
responds to the zero-point displacement of free lattice os- 
cillators. We use Wo as the unit of energy. 

In the same way as for polarons [T], two fundamen- 
tally different kinds of correlative behavior between the 
electron and lattice subsystems can be distinguished for 
HH bipolarons. Adiabatic correlations are described by an 
electron pair that instantaneously adjusts to the motion 
of the lattice deformation field. On the other hand, during 
nonadiabatic processes, the electrons temporarily detach 
from the lattice field. 

It is frequently assumed for [/ = that the parameter 
t/uJo is sufficient to distinguish between the regime domi- 
nated by the nonadiabatic correlations from the one dom- 
inated by the adiabatic correlations. However, it is impor- 
tant to stress that this distinction should also account for 
the strength of the electron-phonon coupling g/uJo- That 
is, the adiabatic correlations develop only for sufficiently 
strong couplings, whereas the weak-coupling limit, irre- 
spectively of t/ujQ, involves purely nonadiabatic dynamics 
[1]. For a fixed t/uJo, the latter can always be reached by 
decreasing g/uio- 

One important ingredient of the bipolaron physics are 
lattice coarsening (discreteness) effects. The role of these 
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effects is determined by the electron-phonon correlation 
length df)p. By analogy with the polaron case [1], the dis- 
creteness of the lattice deformation pins small bipolarons 
{dbp ~ a), whereas in the opposite limit of large bipolarons 
{dbp ^ a), the continuum approximation can be invoked. 

3 Adiabatic approximation 

In the adiabatic limit the electron part of the bipolaron 
wave function behaves as if it commutes with the lat- 
tice kinetic energy ^ , depending parametrically on time, 
through the lattice deformation, 

h(u))=^r;„,„cl_^4jO), (3) 

n , m 

where, for the sake of brevity, the lattice deformation is 
denoted by an TV-dimensional vector u = {un}, with N 
being the number of lattice sites {N — >■ oo). 

For the singlet spin configuration, the two-electron wave 
function ^ satisfies 77„,rn = ilm,m whereas for the triplet it 
satisfies rjn^m = —ilm,n- Furthermore, in the [/ — ^ oo limit 
the singlet solution becomes degenerate with the triplet 
solution. Since in the current, as well as in other works 
|13I20| no indication of stable U = oo singlet bipolarons 
is found, it is expected that the triplet bipolarons are un- 
stable in the whole parameter space of the ID HH model. 

With appropriate normalization m^n mVn,m = 1, 
the expectation value of the Hamiltonian ([Ij with respect 
to |''7(u)) is obtained as 

Had =f + ujQU^+eAD{n), (4) 

where the adiabatic electron energy EAoi^) is given by 
the ground-state (i = 0) solution e(''^(u) of 

~t {Vn+l.rn + (5) 

The last two terms in equation (j4|) are functions of u. This 
means that they can be interpreted as the lattice potential 
energy, henceforth referred to as the adiabatic potential 
Uad (u) . Sad (u) defines the change of the lattice potential 
energy due to the adiabatic correlations, with respect to 
the free-lattice case. The excited-state solutions (i > 0) of 
equation ([S]) may be used to analyze nonadiabatic effects. 

3.1 Adiabatic Holstein {U = 0) bipolaron vs. polaron 

For [/ = 0, the electron part of the adiabatic wave func- 
tion can be factorized as a product of single-electron wave 
functions. Thus, the energy of the electron subsystem is 
given by the sum of single-electron energies e),/ , the latter 
being solutions of 

^el Vn = ~t {r/n+1 + Vn-l) " 2^ Vn , (6) 

with r]n the single-electron wave function. 



At this point, it is convenient to make use [T] of the 
sum rule 

Un = Nel —Un , M„ = 1 , (7) 

'^O n 

where N^i is the number of electrons in the system. This 
sum rule follows from the fact that the homogenous q = 
lattice mode couples only to the total electron density [53] , 

which is fixed. Introducing {t„ and A'^^^ = e"^i iN^iSp in 
equation ^ as rescaled quantities, 

^•■'^ rjn = -TT— r {ln+1 + Vn-l) ~ 2 Un Tjn , (8) 
JVeZ A 

with parameters Sp and A given by £p = g'^ /ujq and A = 
£p/t, the adiabatic potential UAoi'^) rewritten in terms 
of those quantities takes the form 

UADiu)=ep (N^iY,ul + NeiJ2^^'Hun,N,iX) \ . 

\ n i=0 / 

(9) 

Here, the summation over i involves the lowest N^i occu- 
pied single-electron states of equation (|S]). 

Equation ([9]) is general and valid for any N^i . In partic- 
ular, for the bipolaron {Nei = 2) and the polaron {Nei — 1) 
case, only the ground (i = 0) state of the electron spec- 
trum ([8]) contributes to Uad{u). It is singly occupied for 
the polaron and, due to the spin degeneracy, doubly so for 
the bipolaron. An important consequence of this property 
is that Uad{u) for N^i = 2 exhibits the same behavior as 
for Nel — 1 with four times larger Sp and twice as large A, 

£p <^ 4ep , A ^ 2A , g ^2g , t ^ 2t . (10) 

Thus, the adiabatic Holstein bipolaron problem can be 
simply mapped through ([TU|) to the adiabatic Holstein po- 
laron problem. This useful property has not been noted 
previously. 

The quantity in the brackets on the right hand side 
of equation ^ depends only on one parameter, A. As 
known from the polaron theory |31], A defines the adia- 
batic electron-phonon correlation length (the polaron size) , 
dpoi/a w 1 -I- 2/A. Using equation PT?|) . one finds that for 
the same parameters the adiabatic Holstein bipolaron is 
always smaller than the polaron, i.e., dbp /a ~ 1 + 1/A. This 
also means that the pinning effects due to the discreteness 
of the lattice field are stronger for bipolarons than for po- 
larons. In addition, the bipolarons are characterized by a 
doubled lattice deformation in equation ([7]) , which renders 
them heavier than polarons. 

3.2 Nonadiabatic corrections for (7 = 

While the simple mapping (|10p exists between adiabatic 
bipolarons and polarons, such an analogy is absent for 
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nonadiabatic corrections involving the excited states of the 
adiabatic electron spectrum ([5|). For polarons, the ground 
(i = 0) and excited states (i > 0) of the adiabatic electron 
spectrum are given by 



(11) 



whereas the bipolaron case involves one- and two-electron 
excitations, 



2£. 



yl«(u„,2A)+^(^')(u„,2A)' 



(12) 



The comparison of the two spectra in equations 
and (jl2p reveals an important physical property: for the 
same parameters the bipolaron is "more adiabatic" than 
the polaron, because the former is characterized by a larger 
gap in the electron spectrum between the ground and 
excited states. In this respect, it is instructive to consider 
two opposite limits, the large and the small (bi)polaron 
limit, corresponding to A <C 1 and A ^ 1, respectively. 
For large adiabatic polarons, d{,p 3> a, the gap in the elec- 
tron spectrum can be evaluated in the continuum ap- 
proximation |35ll| . which gives = Ep A/4. According 
to equations PT|) and (H^), for the same parameters a 
four times larger gap is obtained in the bipolaron case, 
^?)^ = Ep A. In the small (bi)polaron limit, dtp « a, the 
gap in the electron spectrum is independent of A, being 
two times larger for bipolarons {A'^p = 2Z\55°' = 4ep). For 
arbitrary A, it can easily be verified numerically that the 
ratio /AP°^ lies between the two limiting behaviors dis- 
cussed here, 2 < Z\^p/Z\p°' < 4. 



3.3 Adiabatic bipolarons for [/ 7^ 

In the adiabatic regime the main contribution to the bipo- 
laron binding energy is described by the adiabatic po- 
tential Uad{u), while the kinetic part of the energy, as 
well as nonadiabatic corrections, contributes much less to 
the total bipolaron binding energy. For this reason, some 
important properties of the adiabatic bipolarons may be 
obtained simply by calculating the equilibrium lattice de- 
formation, corresponding to the minima of ?7ad(u). The 
approach that yields the bipolaron behavior from these 
minima is hereafter referred to as the static adiabatic ap- 
proximation (SADA), reflecting the fact that only the adi- 
abatic equilibrium point in the lattice configuration space 
is being considered, while the dynamics are neglected. 



3.3.1 Small bipolarons 

The effects of the Hubbard repulsion on the formation 
of the small bipolarons have been discussed in numerous 
works. For U = and A large, the pair of electrons and the 
accompanying lattice deformation localize to a single lat- 
tice site, forming a so called SO bipolaron |4I13I16| . With 
the weak Hubbard repulsion switched on, the condensa- 



tion energy is given by A^^ 



2er, — U. However, for a 




Figure 1. (Color online) The SADA condensation energy of 
the small (A = 2) and large (A — 1/2) bipolaron as a function 
of U/sp. For the small bipolaron case, different symbols are 
used for SO, 5*1 and 5*2 condensation energies. 



the same lattice site. The resulting bipolaron consisting 
of two weakly overlapping small polarons at the nearest 
neighbor sites is usually referred to as an 5*1 bipolaron 
By keeping only the leading corrections in 1/C/ and 1/A 
small, the SADA condensation energy of the 5*1 bipolarons 
takes the form 



aSI 



4tyu - Sp/X^ 



(13) 



stronger repulsion U > 2ep the electrons tend to avoid 



In fact, A^p can easily be interpreted by starting from 
the U = 00 limit. For U = 00, placing two small polarons 
next to each other divides the energy gain associated with 
the adiabatic spreading of the small polaron by half. This 
repulsive effect is described by the 1/ A^ correction to the 
small-bipolaron binding energy Sp, i.e., by the second term 
in equation Returning to the finite U case, the SI 

bipolaron is stabilized for U < Asp by the superexchange 
energy, given by the first term in equation ([T^ . 

In Figure [T] the binding energy of the small bipolaron 
(A = 2) is shown as a function of U/Sp by curves with 
symbols. The curves are normalized by the U = value 
of Abp. For A = 2, the SADA transition in Figure [T] be- 
tween the SO (circles) and SI (squares) bipolarons occurs 
for U/Sp w 1.8. The second SADA transition in Figure [T] 
takes place for U/Sp ~ 3.8, involving a transition between 
51 (squares) and 5*2 (diamonds) bipolarons. 5*2 denotes a 
bipolaron consisting of two weakly bound small polarons 
at next-nearest neighbor sites. Increasing U further, the 
SADA gives weakly overlapping small polarons at increas- 
ing distances (Si, i > 2), with a vanishing binding energy. 
By calculating the minima of the adiabatic potential as 
a function of U to the leading correction in 1/A small, it 
may be shown |38j that two polarons become unbound for 
a critical value of the Hubbard repulsion Uc- In particu- 
lar, for 1/A = this value is given by Uc = 12 £p, with 
Uc shifting towards larger values as A decreases away from 
the atomic limit 1 / A = |3H] ■ 

Because of strong lattice coarsening effects, the 50 and 
51 bipolarons, corresponding to two minima (local and 
global) of the adiabatic potential Uad{u), are separated 
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by a substantial energy barrier. Consequently, at the tran- 
sition between different Si bipolarons {i ~ 0,1,...) the 
SADA condensation energy Ai,p in Figure [1] exhibits a 
pronounced singularity in its slope. This behavior should 
be contrasted to the large bipolaron case dbp ^ a dis- 
cussed in Section 13.3.21 for which the lattice coarsening 
effects are negligible and the SADA condensation energy 
Abp behaves smoothly [see the A = 1/2 case in Figure [1]. 

In general, kinetic and nonadiabatic contributions are 
expected to smear any singular behavior of the exact con- 
densation energy Abp. Indeed, unlike in Figure [1] Ai,p cal- 
culated numerically by the RCSM in Section H] always ex- 
hibits a smooth crossover between different Si bipolarons. 

3.3.2 Large bipolarons 

The problem of the large HH bipolaron has attracted much 
less attention in the literature than has the problem of the 
small bipolaron. In the continuum approximation, appro- 
priate for dbp ^ a, Uc ^ 2.5 Sp has been reported |36l37j 
as the critical strength of the Hubbard repulsion above 
which the adiabatic bipolaron is unstable with respect to 
forming two separate adiabatic polarons. This value has 
been obtained by a variational technique involving a prod- 
uct of single-electron wave functions. 

Here, in the context of the SADA, the condensation 
energy of the large adiabatic HH bipolaron is studied by 
calculating the exact minima of Uad{'^)-, without assum- 
ing in advance any particular functional form of the adia- 
batic electron wave function. 

While the behavior of the bipolaron condensation en- 
ergy Abp in the small bipolaron limit 1/ A — >■ may be 
approximately obtained |38j from the leading corrections 
in 1/A small, the general case encompasses the summation 
of the whole 1/A, 1/U expansion of the adiabatic binding 
energy. Within the SADA, this summation is performed 
numerically, by using an iterative scheme proposed in ref- 
erence [4], supplemented with appropriate modifications 
necessary to preserve numerical stability in the large bipo- 
laron limit. Following this procedure, the stability of large 
adiabatic bipolarons is established for much larger values 
than predicted before. 

In Figure [1] the condensation energy Abp of the large 
adiabatic polaron is shown as a function of U/Sp (full 
curve). Comparing A{U)/A{Q) for the small and large 
bipolaron cases (A = 2 vs. A = 1/2), one observes that 
the condensation energy of the large adiabatic bipolarons 
drops more slowly. This also means that the critical value 
of the Hubbard repulsion Uc for the bipolaron stability 
increases as A decreases, shifting Uc to much larger values 
[Uc > 12ep). However, as it may be seen from Figure [1] 
irrespectively of A, which defines the size of the bipolaron, 
Abp is rapidly suppressed by U/Sp, becoming exponen- 
tially small for U/ep> 4. 

Further insights into the formation of large bipolarons 
as a function of the Hubbard repulsion U can be obtained 
by examining the SADA lattice deformation. In Figure [5] 
this deformation is plotted for A = 1 /2 and various values 
of U/2ep = i, < i < A. With increasing U, the large 




Figure 2. The SADA lattice deformation of large bipolarons 
(A = 1/2) as a function of U, showing the dissociation of the 
bipolaron in the adiabatic limit. 

bipolaron in Figure [2] progressively transforms into a pair 
of weakly bound large polarons at increasing distance. As 
may be seen from Figure [3J the lattice deformation at the 
central bipolaron site n = is the largest for U/ep < 4. On 
the other hand, for U jsp > 4 the largest deformation is 
found for non central sites n 7^ 0. A close inspection shows 
that, for some values of U , the center of the symmetry 
of the SADA lattice deformation is in fact between two 
lattice sites, e.g., between the n = and n ~ 1 sites 
for U/2ep = 3 and U/2ep = 4 in Figure [H The same 
symmetry of the lattice deformation is found in the small 
bipolaron limit for the 5*1 bipolarons. 

When Abp obtained by the SADA acquires small val- 
ues, corrections associated with the lattice kinetic energy 
are decisive for the stability of the adiabatic bipolaron. 
The results start to be even more intriguing when nona- 
diabatic contributions assume an important role. In this 
context, the interesting case appears when one starts in 
the U ^ 00 limit with two nonadiabatic polarons. By de- 
creasing U the two polarons condense into a bipolaron, 
which for U ~ may be dominated by the adiabatic dy- 
namics. In order to describe accurately such a kind of 
crossovers involving a mixture of adiabatic and nonadi- 
abatic dynamics one has to rely on numerical approaches 
like the RCSM. 



4 Numerical results 

In Section [3] the adiabatic limit was discussed in terms 
of the SADA results. Beyond this, the application of the 
RCSM allows a dynamical quantum description of bipo- 
larons as well as the extension of the current study to the 
whole parameter space. 

It is instructive to start the numerical analysis with a 
comparison to other methods, when the latter are appli- 
cable, so as to establish the accuracy of the RCSM. For 
ujQ = t = g the value of the RCSM ground state (zero mo- 
mentum K = 0) energy is Ebp = —5.420 ujq, which is close 
to the practically exact value of Ebp ~ —5.424 ojq [T3]. For 
the same parameters, high accuracy {Ebp — —5.419 wo) is 
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also achieved by a variational method described in refer- 
ence [6]. 

4.1 Nonadiabatic contributions 

The differences between the nonadiabatic dynamics of bipo- 
larons and polarons become evident by comparing the 
spectra ioi U = through the scaling given by equa- 
tion ([TU)) . For this purpose Figure [3] is used, with RCSM 
curves obtained by varying A ~ Xbp = '2Xpoh while t = 
hp = tpoi/2 is kept fixed. For the right panel of Figure [3] 
the ratio t/uiQ is chosen to be much larger than for the 
left panel in order to contrast the behavior obtained close 
to the adiabatic limit with the regime where the nonadia- 
batic effects play a significant role. For the different choice 
of A scales, both panels exhibit similar band-narrowing ef- 
fects. However, due to the very different values of t/uiQ, 
two different physical mechanisms are involved. 

In Figure [31 the lowest bipolaron band is represented 
by the gray area, with boundaries defined by the K = 
and K ~ n states. This bipolaron band is compared, using 
the scaling in equation ([TU|) . to the energy of the lowest 
polaron band, whose boundaries are given by the full thick 
K = and K = n curves. The dashed curve in Figure |3] is 
twice the energy of the polaron ground state, plotted for 
A = Xbp and t = ti,p. It defines the threshold energy for 
the stability of the bipolaron K = (ground) state. 

For Xpoi = 2Xbp and tpoi = 2tbp, in the absence of 
nonadiabatic contributions, the bipolaron and the polaron 
bands in Figure [3] should be the same. However, the left 
panel {tpoi — 2tbp — coq) in Figure [3] clearly shows a larger 
bandwidth of the polaron band, indicating that the nona- 
diabatic contributions are more efficient in delocalizing 
the polarons. As the decrease of the adiabatic gap in the 
electron spectrum Z\,, leads to an increasing importance of 
nonadiabatic effects, the differences between the bipolaron 
and the polaron band become more pronounced towards 
the weak-coupling (left) side of the left panel in Figure [3] 
(Ajj < Wo). The explanation is quite simple. In nonadia- 
batic processes the electron is temporally detached from 
the lattice deformation. Since two electrons have to work 
cooperatively in order to nonadiabatically delocalize the 
bipolaron. it is not surprising that the polaron delocalizes 
nonadiabatically more efficiently through single electron 
processes. 

On the contrary, for the large t/ujQ used in the right 
panel of Figure [3] (ipoi = 2tbp = GAujq), the differences be- 
tween dispersions of the polaron and the bipolaron bands 
are hardly seen. In other words, the dynamics is almost 
completely adiabatic. The significant band-narrowing in 
the right panel of Figure [3] is governed by lattice coars- 
ening effects that become stronger as A increases, rather 
than by the change in the nature of the electron-phonon 
correlations (i.e., adiabatic vs. nonadiabatic), as it is in 
the left panel of Figure [31 

The qualitative difference between the two (bi)polaron 
band-narrowing mechanisms in the two panels of Figure [31 
may be argued further from the behavior of the bipolaron 
binding energy Ebp. This energy in Figure [31 corresponds 




Figure 3. (Color online) The RCSM lowest band, bounded by 
K = and K = ir states, for the Holstein bipolaron (gray area) 
and polaron (full thick curves) are compared as a function of 
Xpoi = 2Xbp for tpoi = 2tbp = (jJo (left panel) and 64 ujo (right 
panel) fixed (note different A and energy scales in two pan- 
els). The dashed curves are twice the RCSM polaron ground 
state energy for A = Xbp, defining the energy threshold for the 
bipolaron stability, (cjq is used as the unit of energy.) 

to the energy difference between the minimal energy of 
two free electrons —Atbp and the bottom of the bipolaron 
bands. As one may observe, Ebp takes very different values 
in the left and right panels of Figure [31 In particular, the 
small (large) binding energy in the left (right) panel of 
Figure [3l directly indicates the small (large) gap in the 
adiabatic electron spectrum oc Ebp, discussed already 
in connection with equations (fTTj) and (jl2p . For < ojq 
nonadiabatic dynamics prevails, while Z\,, > wq represents 
the opposite, dominantly adiabatic behavior. 



4.2 Bipolaron band structure 

Depending on parameters, excited coherent bipolaroir bands 
may emerge below the phonon threshold for the incoher- 
ent scattering. In fact, due to the energy constraint in sit- 
uations when the low-frequency coherent bands are suffi- 
ciently narrow, coherent bands can also be found above the 
threshold energy. Namely, with narrow bands at the bot- 
tom of the spectrum, some parts of the spectrum above the 
phonon threshold may remain inaccessible to incoherent 
phonon excitations that add wq (optical phonon energy) 
to the total energy of the system. Analogous behaviors of 
the polaron spectrum have been found in investigations 
reported in references |40I41| . 

In a manner similar to that discussed for the lowest 
band in connection with Figure [31 under the substitution 
of parameters (|TU|) any differences between the excited po- 
laron and bipolaron bands in the [/ = case should be 
attributed to nonadiabatic effects. Since these differences 
do not bring any essentially new behavior and since the po- 
laron band structure as a function of the coupling constant 
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g has been extensively reported upon previously |1I27I32| . 
we turn instead to the role of the Hubbard repulsion. 



4.2.1 Relation to soft normal modes 

By revealing excited bipolaron bands, the RCSM is able to 
provide a detailed perspective of various aspects of bipo- 
laron formation. The excited bands are associated with 
the adiabatically softened phonon modes that move along 
the lattice with the bipolaron. Thus, by analogy with the 
polarons [27] , the excited bipolaron bands serve as finger- 
prints of adiabatic correlations. Generally speaking, the 
strongest adiabatic correlations should be expected in the 
low-frequency part of the bipolaron spectrum, for which 
the corresponding lattice deformation is the slowest. The 
absence of the excited bipolaron bands below the phonon 
threshold therefore indicates that, for all frequencies, the 
dynamics is nonadiabatic. Namely, the electrons detach 
nonadiabatically from the phonon cloud too frequently 
and the fiuctuations of the lattice at different sites remain 
adiabatically uncorrelated (non-softened) by electrons. In- 
deed, in the weak-coupling regime (Z\^ < ujq) one always 
finds only the lowest coherent bipolaron band below the 
phonon threshold. 

Due to the softening of the phonon modes, the bipo- 
laron band structure starts to be particularly intriguing 
upon approaching the adiabatic limit. More specifically, 
by varying U, the bipolaron band structure transforms 
from the polaron-like behavior at J7 = [in the sense of 
equation PH)) ] to the large U behavior, when the bipo- 
laron consists of two weakly overlapping polarons sharing 
the same center of mass coordinate. This gradual trans- 
formation of the band structure is examined in detail in 
Figures I2E1 

The four panels in Figure 2] show the SADA (equi- 
librium) lattice deformation u„ and the lowest four adia- 
batic (soft) normal bipolaron modes. The first two panels 
in Figure [4] correspond to the entry and exit values of the 
parameter U, as a function of which the bipolaron band 
structure is plotted in Figure [5] The two remaining panels 
in Figure 2] correspond to the entry and exit values of the 
parameter U for which the band structure is shown Fig- 
ure [5] and discussed in Section 14.2.31 in the context of the 
large U limit. 

As in Figure [U the dissociation of the bipolaron as a 
function of U can easily be tracked in Figure U from u„. 
The parameter A in Figure 3] is approximately twice as 
large as in Figure [2] (A = 0.98 vs. A = 1/2). Consequently, 
the bipolarons in Figure |4] exhibit significant lattice coars- 
ening effects, being set by this particular choice of A in the 
vicinity of the crossover regime between the large and the 
small adiabatic bipolarons. 

The normal adiabatic modes of the bipolaron lattice 
deformation field, shown in Figure 01 are obtained using 
the harmonic approximation for the adiabatic potential 
Uad in equation (jH]). They may be distinguished accord- 
ing to the number of nodes and their parity. Namely, the 
modes are even or odd under reflection with respect to 
the center of the equilibrium lattice deformation ed n ~ 0. 
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Figure 4. (Color online) The SADA lattice deformation and 
the lowest four (soft) normal modes as a function of U (A = 
0.98, t/ujo = 200). 



Depending on the parameters used, the energies of the 
modes with different parities may cross. 

Although calculated by breaking the translation sym- 
metry of the adiabatic bipolaron problem, the normal modes 
in Figure 3] give a clear meaning to the internal structure 
of the low-frequency adiabatic dynamics. That is, by con- 
sidering the relationship between the normal mode coor- 
dinates and the adiabatic coordinate that corresponds to 
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the translational bipolaron motion along the lattice, one 
may distinguish between the local dynamics orthogonal 
to the motion of the center of mass and the translational 
motion of the center of mass itself. 



4.2.2 Band structure for small U 

In the small U regime, the main effect of the Hubbard re- 
pulsion is to increase the size of the bipolaron. For this 
reason the bipolaron band structure in Figure [5] bears 
many resemblances to the small to large polaron crossover 
|1I32| as A is varied. Starting with narrow bands on the left 
(small bipolaron) side of Figiu-e [5l as [/ is increased the 
widths and the distances between various excited bands 
become comparable. As U is increased further the band 
structure on the right side of Figure [5] develops the large 
bipolaron behavior. 

The bipolaron bands in Figure [5] are given in terms of 
9 states with different momenta K, K ~ m x 7r/8, with 
< m < 8. All the energies are shifted by the ground state 
energy of the bipolaron. In order to simplify the analysis, 
additional bands below the phonon threshold associated 
with higher normal modes are not considered in Figure [S] 
(e.g., iV3, N4: modes shown in Figures and |4}d, and 
higher modes). 

The detailed understanding of the band structure in 
Figure [5] involves a distinction between two basic effects 
|1I32| . The first explains the bandwidths, and is related 
to the so called Peierls-Nabarro (FN) barrier Z\pjv- The 
second explains the hybridization between excited bands 
and is related to the effective coupling between normal 
modes. 

Regarding the role of the FN barrier, shifting the large 
adiabatic bipolaron across the unit cell from its exact equi- 
librium position has a very small energetic cost ApN ^ 
loq. Namely, because the lattice coarsening effects are sup- 
pressed for dbp 3> a, the shape of the large bipolaron is 
almost preserved at any point of the minimal energy path 
for the bipolaron translation that connects the exact min- 
ima of the adiabatic potential Uad{u). This effect may be 
observed in the frequency of the pinning P mode, which 
vanishes in the Zip^r — )■ limit. Consequently, the large 
bipolaron moves along the lattice almost as a free parti- 
cle. On the other hand, due to the lattice discreteness, the 
shape of the small bipolaron changes substantially within 
the unit cell along the minimal energy path for the trans- 
lation. Therefore, for small bipolarons the FN barrier is 
large, Zip^r S> ujq, which results in the very narrow bipo- 
laron bands in the spectrum observed in the left part of 
Figure \5\ Their positions in the spectrum correspond to 
the excitation energies of the normal modes shown in Fig- 
ure m At the left side of Figure [S] the bands are denoted 
accordingly, e.g., the band associated to the simultaneous 
excitation of the pinning and breather mode is denoted by 
BP. 

The second effect that requires consideration is the ef- 
fective coupling between the normal modes. The lowest 
even (breather) mode, denoted by B in Figures 0^ and 
HJj, involves vibrations of the bipolaron size, whereas the 




1 2 u/8 3 

P 

Figure 5. (Color online) The bipolaron band structure as a 
function of U {t/ujo = 200, g = 14, A = 0.98). The lowest and 
excited bands corresponding to the excitations of the pinning 
P and breather B normal modes are shown. 

lowest odd (pinning) mode, denoted by P, displaces the 
center of mass. The change in the shape of the small bipo- 
laron along the minimal energy path results in a strong 
effective coupling between the P and the B mode |ll32j . 
This explains the strong hybridization of the excited bands 
in the central part of Figure [5j 

For the large bipolarons in the right part of Figure [SJ 
gaps between various bands (associated with different ki- 
netic energies of the translational motion) close as Appf 
0, whereas the hybridization between excited bands in- 
volving different degrees of freedom (e.g., P and B) be- 
comes weak. B on the right side of Figure [S] denotes the 
K ~ state at the bottom of the band associated with 
the excitation of the breather mode. The position of this 
A' = state in the spectrum is approximately given (up to 
the kinematic effects [33]) by the frequency of the breather 
B mode of the large bipolaron, shown in Figure |3}d. 



4.2.3 Band structure for large U 

Turning now to large U, it is convenient to discuss the 
bipolaron properties in terms of two overlapping polarons. 
Starting with Figures 21 instead of B for breather, P 
for pinning, etc., it is useful to introduce a new cate- 
gorization of the normal modes, iVp, N^, Nb, Na- Here, 
Nji and Nr can be understood as bonding and antibond- 
ing combinations of the two pinning modes correspond- 
ing to the two overlapping polarons. The normal mode 
Nji is odd and displaces the center of the bipolaron mass, 
whereas the mode Nr is even and describes vibrations of 
the distance between two polarons. Similarly, the next two 
normal modes Nf, and Na may be understood as bond- 
ing (even) and antibonding (odd) combinations of two 
breather modes, representing the in phase and antiphase 
vibrations of the polarons' sizes. 
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Figure 6. (Color online) The large bipolaron band structure 
for large U {t/uJo = 200, g = 14, A = 0.98, as in Figure [5]). 
Beside the bands associated to the translation motion of the 
bipolaron with zero point motion of the normal modes, the 
bands with excited normal modes Nr and A^;, are shown as 
well. 

The bipolaron band structure for large U is shown 
in Figure [51 Steps in gray shading represent the increase 
in the number of overlapping bands in the spectrum. As 
in Figure [5l all energies are shifted by the ground state 
energy. The band boundaries correspond to the K ~ 
(full curves) and K = tt (dashed curves) states. Which 
of the two K states actually defines the bottom and the 
top of the corresponding band depends on the parity of 
the states. For example, the ground K = and the lowest 
K = TT state in FigureOare even, and define the lower and 
upper boundaries of the lowest band, respectively. A small 
gap (white area) separates this lowest band from the first 
excited band that starts with the K = it state of odd par- 
ity. In the absence of lattice discreteness effects, the gap 
between the lowest two bands closes, and the dispersion of 
the bipolaron states is given simply by Ei,p{K) (x /rribp, 
where rubp denotes the bipolaron effective mass. The states 
associated with the momenta K — it and X = 27r (corre- 
sponding to the Nb, coordinate) are represented separately 
in Figure [HI 

Unlike in Figure [51 where the effective mass m^p of 
the bipolaron decreases on increasing the bipolaron size 
with [/, in Figure[S|C/ has the opposite effect. Namely, the 
bandwidth of the lowest band clearly decreases with [/, 
indicating that the polaron pair becomes heavier as the 
mutual distance between the polarons increases. 

Beside the bands associated with the increasing kinetic 
energy of the joint motion of two polarons along the lat- 
tice, in Figure [Hi additional bands associated to two even 
normal modes Nj. and Nf, are shown, with the K = Q 
states denoted according to the nature of the excitation 
involved. In particular, it may be seen from the K = Q 
state denoted by Nr in Figure [Hi that the frequency of the 
Nr mode decreases with U . This behavior is expected on 
the basis that the restoring force for the vibrations of the 



distance between two polarons vanishes when the bipo- 
laron dissolves into two unbound polarons. 

It is also worth noting that, by sharing the same parity, 
the NrNr state (double excitation of the Nr mode, K = 
0) and the Nj, state (single excitation of the Nt mode, 
K = 0) , anticross in the central part of Figure [5] That 
is, the K = state denoted by Ni, on the left side of 
Figure [5| changes its nature, being dominantly a double 
excitation of the Nr mode on the right side of Figure [5] In 
other words, the hybridization between bands associated 
to the Nf, and Nr excitations of the bipolaron occurs in 
Figure [S] However, in contrast to the hybridization due 
to lattice coarsening effects in Figure [SJ in Figure [Sj one 
observes only a weak effect of the kinematic j35] origin, 
without abrupt changes in the dispersion properties. 



4.3 Light bipolarons with significant condensation 
energies 

While the adiabatic limit ^ uiq involves large lattice 
deformations that make the bipolaron heavy, the weak- 
coupling limit corresponds to the opposite situation. In 
this respect, it is interesting to consider which values of 
U give the most stable, light bipolaron solutions. The ex- 
istence of light bipolarons with significant condensation 
energies has attracted particular attention in the context 
of bipolaron superconductivity |42| and findings that in- 
dicate the importance of the electron-phonon interaction 
in higli-Tc materials |43| . 

In the HH model, one finds [131316) that light bipo- 
larons with significant binding energies exist when the rel- 
evant energy scales governing the bipolaron dynamics are 
comparable, t ~ Z\,; ^ wq. It is emphasized here that this 
specific regime of parameters corresponds to the crossover 
between the weak-coupling limit and the regime of small 
bipolarons characterized by strong adiabatic correlations 
and lattice coarsening effects. In other words, light bipo- 
larons with significant binding energies involve a subtle 
balance of parameters. 

Figure [7^ shows the RCSM bipolaron effective mass 
rubp for three values of the Hubbard repulsion, U/ujo ~ 

0, 2,4. The polaron effective mass nipoi (dashed curve) is 
plotted for comparison. A weak renormalization of rupoi 
indicates the weak-coupling regime for the polaron case, 

1. e., the dynamics are fully nonadiabatic. On the other 
hand, in contrast to polarons, the bipolaron spectrum (not 
shown) exhibits excited bands below the phonon threshold 
for the set of parameters investigated in Figure [71 meaning 
that the adiabatic correlations are significant. 

The RCSM bipolaron ground-state energy (full curves) 
is shown in Figure [TId for the same set of parameters used 
in Figure [7^. In the small (5*0) bipolaron limit {SpX ^ 
Ep ^ uJo and Sp > U), corresponding to the right side of 
Figure [Tb, the ground-state energy approaches asymptot- 
ically the values plotted by the thin dashed curves and is 
given by 

Ebp « -4ep - Ep/X'' + U . (14) 



10 



O. S. Barisic, S. Barisic: Bipolarons and polarons in the Holstein-Hubbard model: Analogies and differences 



0.5 



0.25 



1 


'1 

a) 






\U=4 






— ..^ ^ — 1 





1.4 1.6 1.8 2 




1.4 1.6 1.8 2 



Figure 7. (Color online) a) The RCSM bipolaron effective 
mass for U /ujq — 0, 2, 4. b) The RCSM bipolaron ground-state 
energy for the same parameters. 



The first term is the energy of two electrons localized at 
the same site (5*0 bipolaron), while the second term is 
the energy gain due to the adiabatic spreading of the 50 
bipolaron to the two neighboring sites. Upon substitution 
of parameters (fTO)) . the first two terms in equation (jl4p 
represent the energy of the small adiabatic polaron cal- 
culated to the leading order in small 1/A. The effects of 
small Hubbard repulsion in equation (jl4p are taken with- 
out 1/A corrections, as if two electrons were permanently 
sharing the same lattice site. 

In the crossover towards weak couplings (central part 
of Figure [TId) , significant deviations from equation (jl4p 
start to occur. Namely, with decreasing g and/or increas- 
ing U, the spreading of the lattice deformation renders 
bipolarons lighter and simultaneously suppresses the gap 
in the adiabatic electron spectrum Z\^. Such suppression 
introduces significant nonadiabatic correlations. 

The range of nonadiabatic correlations increases rad- 
ically in Figure [Tb for states close to the threshold en- 
ergy for the bipolaron stability, when tire bipolaron dis- 
solves nonadiabatically into two polarons. The description 
of such weak and long-ranged nonadiabatic correlations is 
quite approximate within the RCSM (and other applica- 



ble methods) and small inaccuracies become notable for 
^bp 0. In particular, instead of approaching asymptot- 
ically the thick dashed Atp = curve from below in Fig- 
ure [7)3, tire RCSM ground-state energy curve for U — Aluq 
intersects it. Yet, for this particular regime of parameters, 
the expected relation rubp ~ 2mpoJ is obeyed in Figure [TK- 
This result shows that the overall RCSM picture of the 
bipolaron dissociation for Abp — >■ is semi-quantitatively 
correct. 

Table 1. Bipolaron RCSM condensation energy as a function 
of U for two values of the effective mass. All the energies are 
in units of loq {t = 2.5 cjq). 
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[7 = 


0.38 
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0.2 

(7 = 4 
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1.54 


1.66 


1.8 


^bp 


1.12 


0.58 


0.34 



For the regime t ^ Ajj ^ luq, the relationship between 
the condensation energy Ai,p and the effective mass mbp 
is further elucidated in Table [T] One observes that, for a 
given effective mass, the condensation energy monotoni- 
cally decreases with U. This means that the most favorable 
conditions for the formation of very light HH bipolarons 
with large condensation energies are achieved when the 
Hubbard repulsion is negligible. 



5 Entropy considerations 

In the HH model the effective interaction between indi- 
vidual bipolarons is repulsive |13| . Therefore, assuming a 
positive condensation energy Abp > 0, the electrons tend 
to condense in pairs. A naive expectation in such circum- 
stances is that the majority of charge carriers condense 
into bipolarons up to temperatures comparable to Abp. 
However, the entropy of the electron-phonon system de- 
pends on the density of charge carriers and the ratio of 
the concentrations of polarons and bipolarons in the sys- 
tem exhibits a more intricate behavior. Therefore, while 
the single polaron physics may be invoked for the dilute 
Abp < limit, the bipolaron problem Abp > requires ad- 
ditional considerations. That is, in the dilute limit, even a 
small temperature T can drive the system from the bipo- 
laronic ground state (T = 0) to the polaron phase because 
the latter is favored by the gain in the free energy through 
the increase of the entropy S. This fact is frequently over- 
looked and, instead, only the binding energy is used to 
estimate the relative concentrations of polarons and bipo- 
larons. 

In this connection it is instructive to consider the small 
(bi)polaron limit, for which the length of local electron- 
phonon correlations d reduces to just one lattice site, d/a ^ 
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1. In this limit, the analysis is considerably simplified be- 
cause of the vanishing overlaps between bipolarons and 
polarons, while the Pauli exclusion principle prohibits the 
double occupancy of lattice sites by electrons with the 
same spin. If there are N lattice sites, the number of per- 
mutations of placing two electrons involved in small po- 
larons on the lattice is N{N — l)/2, whereas the small 
bipolaron can be placed at N different sites. It follows 
that the gain in entropy of forming two polarons increases 
logarithmically with the system size S ^ ln{N). The ex- 
cited states of the bipolaron do not play a fundamental 
role here, since the number of them is finite, independent 
of N and given approximately by dbpAbp/acoo. 

For a given total charge concentration c per site the 
average concentration per site of small bipolarons Cbp and 
the spin-degenerate concentrations of small polarons c-f-j^ 
which satisfy c = 2cbp + c-|- + c^, can be expressed by [IJ 



2ct(1 



) = l-(l-c)i 



c(2 - c) 

(1-C)2 



jM^p ^ (15) 



with /3 being the inverse temperature. Although derived 
previously [H], some important aspects of equation (jl5p 
remain to be emphasized. In particular, for any finite P in 
the dilute limit c — )• of equation ([T5)) . all the charge is 
assigned to the polarons |45) . 

^ l-(l-c)(l + ce^-^^^) ^ 

« (i_e/3A.,) « c ■ (16) 

With the increase of the total charge concentration c, the 
ratio c^/c rapidly decreases provided that /3Z\f,p is large. 
For example, for (3Ai,p ~ 10 and c = 0.05 the ratio Cbp/c 
is close to 0.97. Yet. for smaller condensation energies like 
PAbp = 5, large relative values of bipolaron concentration 
Cbp/c > 0.8 are obtained for c > 0.2. Thus, the observation 
of bipolarons in the dilute limit is possible only for low 
enough temperatures (3Abp ^ 1. Otherwise, only polarons 
will be observed. 

In general, for Abp > and dtp arbitrary, overlapping 
polarons and bipolarons are simultaneously present in the 
system, which complicates the estimation of their ratio as 
a function of doping and temperature. Nevertheless, the 
free energy gain related to the formation of two polarons 
instead of the bipolaron is large whenever the correlation 
length satisfies dbp/a ^ c^^. 



6 Summary 

The current work provides a thorough examination of the 
low-frequency properties that characterize the formation 
of bipolarons within the ID Holstein-Hubbard model. Par- 
ticular attention is devoted to the analogies and differences 
with respect to polarons, since the properties of these two 
kinds of quasiparticles are expected to govern the behavior 
of the electron-phonon system in the dilute limit. 

For both, the bipolarons and the polarons, the deep 
dichotomy in the electron-phonon correlations is funda- 
mentally the same, i.e., adiabatic vs. nonadiabatic. The 



first interesting observation along these lines is that for 
U — the adiabatic bipolarons and polarons exhibit the 
same spectrum under the simple scaling of parameters, 
derived here in equation ([TU]). It is next argued that, for a 
given set of parameters, bipolarons are always more adia- 
batic than polarons. This allows some easy predictions of 
the bipolaron behavior using the already known polaron 
behavior. 

As a function of U , two basic limiting behaviors may 
be distinguished. For U small, the repulsion between elec- 
trons increases the bipolaron size and, depending of pa- 
rameters, one may observe a small to large bipolaron crossover 
that is very similar to the small to large polaron crossover 
when the polaron size is varied through X = /t loq. For 
U large, the bipolaron may be discussed in the picture 
of two overlapping polarons that move together along the 
lattice. As shown here, detailed aspects of the small and 
large U behavior may be easily understood from the bipo- 
laron band structure. The latter, including the previously 
unreported excited bands, is calculated by the RCSM, a 
method that has previously been successfully applied to 
the polaron problem. As for polarons, the excited bipo- 
laron bands are associated with the adiabatically softened 
phonon modes of the moving lattice field. When the weak- 
coupling regime is achieved, the softening effects are sup- 
pressed completely by the nonadiabatic dynamics and the 
bipolaron spectrum below the phonon threshold involves 
only the lowest band, just as in the polaron case. 

For weak electron-phonon couplings the bare interac- 
tion between two electrons can be approximated by an 
instantaneous effective interaction (the frequency depen- 
dence of the phonon propagator may be neglected). Con- 
sequently, it may be determined that the electron pair 
binds for U < 2ep. In the adiabatic limit, due to retar- 
dation effects associated with the lattice, the bipolaron 
condensation energy Abp remains positive for large values 
of U . However, for U > 4ep, Abp stays small in the physi- 
cally relevant part of the parameter space. This property 
is independent of the bipolaron size, characterizing both 
the small and the large adiabatic bipolarons. 

Particularly interesting is the behavior of small light 
bipolarons, whose condensation energy is comparable to 
the bare phonon energy wq. Under these conditions, one 
finds that a fine balance is achieved between the adiabatic, 
nonadiabatic, and lattice coarsening effects. The role of 
the Hubbard repulsion is to decrease the effective mass 
•nibp and to suppress Abp. In particular, for fixed mbp, 
Abp decreases monotonically with U . This implies that 
the light bipolarons are most strongly bound when the 
Hubbard repulsion is negligible. 

Within the HH model, the effective repulsion between 
individual bipolarons suppresses the phase separation at 
low charge densities. Yet, the conditions under which a 
strong fraction of bipolarons can be observed are severely 
restricted by the temperature and total charge concentra- 
tion. For low temperatures pAbp 1, the relative concen- 
tration of bipolarons and polarons changes in favor of the 
former with increasing total charge concentration. On the 
other hand, for a given finite temperature and a vanishing 
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total concentration, the system is driven into a polaronic 
phase due to the loss in entropy associated with the for- 
mation of bipolarons. 
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